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ABSTRACT 

The  pure  initial  value  problem  for  the  system  of  equations 

vt  “  v*x  +  f(v)  -  w 
wt  =»  e(v  -  yw) 

is  considered.  Here  e  and  y  are  positive  constants,  and 
f(v)  «  v  -  H(v  -  a)  where  H  is  the  Heaviside  step  function  and 
a  e  (0,  V2).  Because  of  the  discontinuity  in  f  one  cannot  expect  the 
solution  of  this  system  to  be  very  smooth.  Sufficient  conditions  on  the 
initial  data  are  given  which  guarantee  the  existence  of  a  classical  solution 
in  R  x  (0,T)  for  some  positive  time  T  • 
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SIGNIFICANCE  AND  EXPLANATION 


The  most  famous  model  for  nerve  conduction  is  due  to  Hodgkin  and 
Huxley.  However,  a  mathematical  analysis  of  their  model  has  proven  very 
difficult.  The  complexity  of  the  Hodgkin  and  Huxley  model  has  led  a  number  of 
other  authors  to  introduce  simpler  models.  In  this  report  we  consider  one 
such  simplification. 

It  has  been  demonstrated  (experimentally)  that  impulses  in  the  nerve  axon 
travel  with  constant  shape  and  velocity.  Mathematically,  this  corresponds  to 
traveling  wave  solutions.  A  number  of  authors  have  proven  that  the  mathe¬ 
matical  equations  considered  here  do  possess  traveling  wave  solutions. 

Another  property  of  impulses  in  the  nerve  axon  is  the  existence  of  a  threshold 
phenomenon.  This  corresponds  to  the  biological  fact  that  a  minimum  stimulus 
is  needed  to  trigger  an  impulse.  Here  we  prove  some  preliminary  results  which 


will  be  used  in  a  later  report  when  it  is  demonstrated  that  the  equations 


under  study  do  indeed  exhibit  a  threshold  phenomenon. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


LOCAL  EXISTENCE  FOR  THE  CAUCHY  PROBLEM  OF  A 


REACTION-DIFFUSION  SYSTEM  WITH  DISCONTINUOUS  NONLINEARITY 
David  Teraan 

1.  INTRODUCTION 

In  this  paper  we  consider  the  pure  initial  value  problem  Cor  the  FltzHugh-Nagumo 
equations 

(1*1)  vt  -  v^  +  f(v)  -  w 

wt  -  e(v  -  YW)  , 

the  initial  data  being  (v(x,0),  w(x,0))  -  (^(x),0)»  Here  e  and  y  are  positive 

constants.  These  equations  were  introduced  as  a  qualitative  model  for  nerve  conduction 

[2,5,7]*  We  follow  McXean  [4]  and  assume  that  f(v)  is  given  by  f(v)  *  v  -  H(v  -  a) 
where  H  is  the  Heaviside  step  function  and  a  6  (0,  V2  )  . 

Note  that  because  f(v)  is  discontinuous  we  cannot  expect  the  solution,  (v,w),  to 
be  very  smooth.  By  a  classical  solution  of  System  (1.1)  we  mean  the  following! 

Definition:  Let  S^,  «  R  x  (0,T)  and  -  {(x,t)  6  ST:  v(x,t)  *  a).  Then 

(v(x,t) ,w(x,t) )  is  said  to  be  a  classical  solution  of  the  Cauchy  problem  (1.1)  in  ST  if 

(a)  (v,w)  along  with  (vx»wx)  are  bounded  continuous  functions  in  S,p, 

(b)  in  Gt  ,  v^,  vfc  and  wfc  are  continuous  functions  which  satisfy  the  system 

of  Equations  (1.1), 

(c)  lim  v( x, t )  «  <^(x)  and  lim  w(x,t)  -  0  for  each  x  6  R  . 

tlO  t+0 

Throughout  this  paper  we  assume  that  *>(x)  »  v(x,0)  satisfies  the  following 
conditions : 
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U)  *<X)  e  c'(«) 

(b)  *(x)  -  v»(-x)  in  R 

(c)  ^(Xq)  •  x  for  some  xQ  >  0 

(d)  V>(x)  >  a  if  And  only  if  |x|  <  Xq 

(•)  (  Xq  )  <  0 

(f)  <?”(x)  is  a  bounded  continuous  function  except  possibly  at 


This  last  condition  is  needed  in  order  to  obtain  sufficient  a  priori  bounds  on  the 
derivatives  of  the  solution  of  System  (1.1). 

In  this  paper  we  prove  that  if  <i(x)  satisfies  (1*2)  then  there  exists  a  classical 
solution  of  the  Cauchy  problem  (1.1)  in  ST  for  some  T  >  0  •  Here  we  give  an  outline  of 
the  proof* 

From  Assumption  (1*2)  we  expect  there  to  exist  a  positive  constant  T  and  a  smooth 
function  s(t),  defined  in  [0,T] ,  such  that  v  >  a  for  |x|  <  s(t)  and  v  <  a  for 
|x|  >  s(t).  Suppose  that  this  is  the  case.  We  then  set  G  -  {(x,t):|x|  <  s(t), 

0  <  t  <  T}  and  let  Xq  be  the  characteristic  step  function  of  the  region  G.  It  follows 

that  if  |x|  ^sCt),  then  (v,w)  is  a  solution  of  the  system  of  equations 


v  «  v  -  v  +  v  -  w 
t  xx  *G 

w  •  e<v  -  -yw) 


<v(x,0),w(x,0))  -  (v? (x)  ,0 ) 


in  sT  , 


Note  that  the  first  equation  in  (1.3)  is  similar  to  a  nonhomogeneous  heat  equation  while 
the  second  is  just  an  ordinary  differential  equation*  Formally ,  the  solution  of  (1.3)  can 


be  written  as: 


v<x,t>  -  J"  K(x-C,t)V»(5)de  +  /*  dT  fB(T)  K(x-e,t-T)d5 

0  -s(  t) 


-  fdT  r  K(x-C,t-T)w(  5»T)de 
o  —00 


w(x,t)  *  £  e  GYt  /tee^rn  v(x,n)dn 
0 


Hare  K(x,t)  -  - - - r? 

2k  ' 2  t  /2 

equation 


-x^/4t 

e  is  the  fundamental  solution  of  the  linear  differential 


“  *xx  '  ♦ 


(1.5) 

Setting  x  *  s(t)  in  <1.4 )  we  find  that,  formally,  a(t)  must  satisfy  the  integral 
equation 

ft^_  fS(T) 


(1.6) 


a  -  /  K(s(t)-e,t)*<C)<H  +  /  dx  J  1  1C(s(t)-C,t-T)de 

-OD  0  -•(  T) 


-  J*  dx  r  K(s(t)-5.t-T)*r(5,T)d5. 
o  — ® 


Using  an  iteration  procedure,  we  prove  the  existence  of  functions  v( x,t) ,w( x,t) , 
and  s(t)  which  satisfy  the  Equations  (1.4)  and  (1.6).  We  then  show  that  (v,w)  is  the 
desired  classical  solution  of  the  Cauchy  problem  ( 1 . 1 )  in  ST  . 

We  now  introduce  some  notation. 

Let  if-(x,t)  -  /*  K(x-£,t)  <*>(  £)<H  •  Note  that  ^(x,t)  is  the  solution  of  the  linear 

>«s 

equation  i 1  -  5 1  with  initial  datum  ^(x,0)  -  ^(x). 

Suppose  that  a(t)  is  a  positive,  continuous  function  defined  in  [0,Tj]  for  some 

>  0  •  Let  z(x,t)  be  a  continuous  function  defined  in  R  x  [0,T  )•  Let 

*<a)(t)  -  j^dx  /a(T)  K(o(t)  -  £,t-x)d£  In  [0,1.1 
0  -a(x) 


and 


r(s)(x,t)  -  /fcdT  f*  K(x-C/t-T)z(  £,  T)d£ 
0  —00 


in  R  x  [0,T^J  • 


Note  that  s(t)  is  a  solution  of  (1.6)  if  and  only  if 

(1.7)  ip(s(t),t)  -  a  -  *(s)(t)  +  T(w)(s(t),t). 

In  Section  2  we  prove  the  properties  of  and  the  operators  ♦  and  r  which  are 
needed  in  the  proof  of  the  local  existence  of  a  classical  solution  of  System  (1.1).  The 


1 


proof  of  local  existence  is  given  in  Section  3 


In  thie  section  we  prove  the  properties  of  ^  and  the  operators  $  and  Y  which  are 

needed  in  the  proof  of  the  local  existence  of  a  classical  solution  of  System  (1*1). 

00  + 

Lemma  2.1.  ijj(x,t)  0  C  (*  x  R  ).  Furthermore,  there  exist  positive  constants 
^t'^2'^3  and  *  such  that  -6^  <  ^(x,t)  <  -6^  and  |^(x,t)|  <  6^  in  the  rectangle 
P  -  (x  -  \,  xQ  +  \)  x  (0,1). 

Proof ;  The  first  assertion  is  a  standard  result  about  solutions  of  Equation  (1.5).  The 
other  assertions  follow  from  the  Assumptions  (1.2e,f).  (See  Friedman  [3],  page  65).// 


Lemma  2.2. 


(2.1) 


Assume  that  a(t)  6  c\o,T).  Then  $(a)(t)  6  c\o,T)  and 
x 

*( a) 1  (t)  -  /  K(a(t)-C,t)dC  +  /  K(a(t)  +  a(T),t-T)ta'(x)  +  a'(t))dT 

-xo  0 


+  J*  K(  a(t)-a(  x)  ,t- x)  t  a1  (  x)-a*  (t)  ]  dx  . 
0 

Proof :  Note  that 

♦  (a) '  (t)  -llm-  [4(a)  (t+e)  -  4(a)(tU 
€+0  € 


-  llm-  [Jt+edx  Ja(T>  (a(t+e)-5,t+e-T)d5 
£+0  6  0  -a(x> 


-  J*  dT  /a(T>  K(a(t)-5/t-T)dC] 
0  -a(x) 


..  1  ff0  .  fa(  t+c)+a(t)-a(  T+e)  r  ^ 

*  lim  -  [/  dT  J  K(  a(  t) t-x)d{; 

e+0  e  e  -a(  x+e)+a(t)-a(t+e) 


+  /fc  dT  J“a(T)  Kla(t)-5,t-T)d5 

0  -a(T+£)+a(t)-a(t+£) 

-  f  dr  /Q<T+E>+a(t)-aU+E)K(a(t)-e,t-T)d5) 

0  -a(x) 


Passing  to  the  limit  we  obtain  (2.1).  // 


Lemma  2.3:  Suppose  that  o(t)  6  C^(0,T)  and  |  o*  (t)|  <  M  in  (0#T)«  Then 

1  ♦(  a) '  (t)  1  <  1  +  Wei^2  in  (0,T)  . 


Proof:  Prom  (2*1)  it  follows  that  for  t  e  (0,T), 


I  ♦(  a) ' (t) I  <  1  +  4M  /*  K(a(t)  -  a(T),t-r)dr 
0 


<  1  +  4M  f  - t I—1 


<  1  +  4MT 


0  2.  2  (t-T)  2 

1/2  •  // 


dT 


Lemma  2.4:  Let  a(t)  be  as  in  the  previous  Lemma.  Suppose  that  for  some 
there  exists  a  constant  such  that 

1a 

|a,(t1)  -  a'(tQ)  I  <  2 

for  each  tg,  t1  e  (p,T).  Then  there  exist  positive  constants  and  K2# 

only  on  p  and  H  ,  such  that 

l*(d)'(t1)  -  +  k2M1T1/2  >'*1  "  t01  ^ 

for  each  t0,  t1  e  (p,T). 

Proof ;  Fix  tg,  t ^  e  (p, T)  •  Then 

x 

♦  (a) 1  (t, )  -  ♦( a) '  (tQ )  -  /°  -  JC(  a(t0 )-5,t0 )] d? 

_X0 

+  £/  V(a(t1)-a(T)»t1-T)(a*(T)-a,(t1)]dT 


K(a(tQ)-a(T)  ttg-t)  (a'  (T)-a'  <tQ)l  d-tj 


[/’ 


K(a(t1  )+a(T)/t^-T)  [a'  ( tl^a1  (t4| ) ]  dr 


-  {  0  K(a(t0)  +  a(T),t0-T)ta'(T)  +  aMt0)ldTj 

-  [A]  +  [B]  +  [Cl . 


e  <o,t) 


which  depend 


-5- 


! 

! 


Sine*  K(x,t)  is  an  infinitely  differentiable  function  of  <x,t)  for  t  >  0,  there 
exists  a  positive  constant  D1  such  that  | [A] |  <  ^  •  Note  that  D1  depends 

only  on  p  and  M  . 

Next  consider  (B]  which  we  rewrite  as 


[B]  -  /  0  Wa(t1)-a(Ttt1-t0),t0-T)<<»,(ml-t0)-a,<t1))dT 


VS 


/  K(a(t  )-a<x),t  -xHa*  (x)-a' (t  ))dx  . 


Note  that, 


aXx+t^)  -  a’^)  -  [a'Crft^tg)  -  a*(x)J 


+  [a'<T)-a'(tQ))  +  [a'  (tQ )  -  a*^)] 


<  la'  ( x)  -  a'(t0)]  ♦  2M1lt1  -  tQ| 


V2  • 


Therefore, 


I  [B]  I  <  [  J°  S(a(tt)-a(t»t1->t9),t0-T)Iaf,(T»t1-t0)-a'(e1)]dT] 

VS 

+  [  /  0IK(a(t1!-a(T+t0-t1),t0-T)  -  Ktcrf^J-aCxJ^-TUta’tTj-a'CtjJldT] 
♦  £  /  °  K(a(t1)-a<T+t1-t())  ,tQ-T)2M1 1 1 tQ I  dtj 


Now, 


[B,]  +  [B,]  +  [B,J 


3J. 


[B  ]  <  2M  J°  - - 17  dT 

VS  2*  2(t<fT)  2 


.  2"  ,tV2_tV2| 

1  i  o  1 

V2 


<02  «S  -  v 


for  some  constant  Dj  which  depends  only  on  M  and  p  •  We  also  have  that 


-6- 


Therefore 


,3  ,*0  ■VV 


l[B251  <  8,1  t  - IX'  U - 17  1 

2  0  2*  ^  (t  -t)  /2 


3  Vo 

<  8M  T  2  |tj  -  tQ| 


<  Ojlt,  -  t0l  /2 

for  some  positive  constant  Dj  which  depends  only  on  M  .  (Since  we  will  eventually 
choose  T  to  be  small  we  assume  throughout  that  T  <  1 . ) 

We  have  shown  that 

I  IB]  I  <  I  tB1J  |  +  |  [Bj]  |  +  |  [Bj]  |  <  D4|t1-t0(  1/2+  2M1T11/2  |t1-tQ|  1/2 
where  depends  only  on  M  and  p  • 

A  similar  computation  shows  that  there  exist  constants  D5  and  Dg,  which  depends 


only  on  p  and  M  ,  such  that 


\  x  n  m  t^2u  |  ^2 


I  tell  <  °5I  VV  2  ♦  D6M1T1  2  • 

In  fact,  this  computation  is  much  easier  since  K(a(t^)  +  a(T)#t^-T)  and 
K(a(tQ)  +  a(T),tQ-T)  are  smooth  functions  of  x  . 

Setting  +  D4  +  and  K2  *  2  +  ,  The  result  follows •  // 

We  now  consider  the  operator  r<z)(x,t)«  In  what  follows  we  assume  that  T  is  some 
positive  constant  and  S,p  *  R  x  (0,T).  We  also  assume  that  a(t),M,p,  and  are  as 

in  the  previous  two  lemmas,  and  set  h(t)  =  H z) ( a(t) ,t) . 


Lemma  2.5:  Assume  that  z(x,t)  e  C1'1(ST)  with  Hz  I  *  Z  .  Then, 

c  * 

i)  h(t)  e  c’fO.T)  , 

ii)  there  exist  a  constant  ,  which  depends  only  on  Z  ,  such  that 

|  h *  < t )  |  <  K3  +  for  t  e  ( 0 ,  T) , 

ill)  there  exist  constants  and  ,  which  depend  only  on  p,  M,  and  Z 

such  that 

IhMt,)  -  h’(t0)|  <  (K4  +  )|t,-  t„| 

for  each  t^,tj  e  ( p,T) . 


-8- 


Proof:  Set  g(x,t)  -  r(z)(x,t).  Then  g(x,t)  is  the  solution  of  the  inhomogeneous 

differential  equation 

-  Sex  "  u  +  z 
u(x,0)  -  0  . 

Since  z  0  C1'1(ST)  it  follows  from  the  Schauder  estimates  (see  [3],  page  65)  that 

g  0  C2'  ^2  (S-)  where  Igl  *  y,  depends  only  on  Z  •  We  set  K  -  Igl  y. 

C2'  *  (B  )  6  C2'  *  (ST) 

Furthermore,  there  exists  a  constant  ,  which  depends  only  on  Z,  such  that 

|gt(x,t)|  <  K3  in  ST.  Note  that  in  ST  , 

|gx(x,t)|  -  |/fc  dt  r  Kx<*-e.t-T)*(t.T)<Sel 
0  -® 

-  |  jSst  J°°  K  (x-£,t-x)z(  £,  r)d£| 

0  -«0  ^ 

-  I  dx  J°°  K(x-C,t-x)z  (C,T)d^| 

0  -«  * 


<  Z  J^dx  f°  K(x-£,t-x)d£ 


<  z 


0  (t-T) 


dx 


<  Z  T 


V 


Now  (i)  follows  because  h(t)  *  g(a(t),t)  where  g  and  a  are  both  smooth 
functions.  (ii)  is  true  because 

(2.2)  hf (t)  -  g  a(t) ,t) a* (t)  +  g  <a(t),t)  in  (0,T) 

x  t 


and,  therefore, 


finally,  it  follow*  fro«  (2.2)  that  for  tQ#t^  «  ( p,T) » 

|h*(tt)  -  h'  (tg)  |  <  |gx(a<t1),t1)  -  gx(  a(t0>  ,tQ>  I  t  <*•  <  | 

+  |gx(a(t0),tQ)  |  la'ttj)  -  a*<t0)|  +  lgt< orftj)  9t^ )  - 

<  K6M|trt0l  +  CT1/S*M1lt1-t0l  K6lt,  -  t0l  V‘ 

(iii)  now  follows  if  ws  set  *4  *  K$(M  ♦  1)  nnd  «  S  .  // 


V^V'V 


-10- 


3*  Local  existence 


We  are  now  ready  to  prove  the  existence  of  a  classical  solution  of  the  Cauchy  problem 
(1.1)  in  ST  for  some  positive  T  •  The  idea  of  the  proof  is  as  follows* 

Let  »0(t)  *  Xq  in  R+  and  suppose  that  for  some  time  Tj  >  0  we  have  defined 
smooth  functions  s^lt)  for  t  ©  1 0 # 3  #  k  *  0,1,***,n.  We  then  let 
( vn(x,t) ,wft(x,t) )  be  the  solution  of  the  integral  equations 

m  *.  »  <T) 

v  (x.t)  -  /  K(x-e,t)v>(5>ae  +  r«*T  /  Ku-e,t-T)d£ 

0  -8  < T) 

n 

(3.1a) 

-  J^dT  f*  K(x-£,t-x)w  (£,r)d£ 

'  n 

0  -• 

(3.1b)  w  (x,t)  *  ee  eYt  /*  eeYT1  v  (x,n)dn  * 

n  o  n 

That  such  a  solution  exists  is  proved  in  Lemma  (3.1)*  We  then  use  the  Implicit  Function 

Theorem  to  define  sn+1(t)  as  the  solution  of  the  equation 

8  ^<t),t>  -*-♦<•  )<t)  ♦  r(w  )  (s  (t)  ,t)  , 
n+ i  n  n  n 

s  (0)  -  x  . 
n  0 

We  show  that  the  sequences  of  functions  {a  (t)},  {v  (x,t) },  and  {w  (x,t)}  converge  to 

n  n  n 

functions  s(t>,  v(x,t),  and  w(x,t).  These  functions  are  shown  to  be  solutions  of  the 
Equations  (1*4)  and  (1.6).  It  is  then  shown  that  (v,w)  is  a  classical  solution  of  the 
Cauchy  problem  (1.1). 

In  what  follows  we  let  sQ(t)  *  x^  in  R+  and  assume  that  smooth  functions  s^tt), 

k  *  0,1,  ♦•*,!*,  have  been  defined  in  [0,T.j]  for  some  >  0  .  Restrictions  on  Tj  will 

be  given  later.  We  assume  that  fi  -  sup  |s*(t)|  <  •  ,  for  k  -  0,1,2,  •••,«•  Por 

teco,^) 

each  p  e  (0,^)  we  assume  that  there  exist  constants  such  that 

I  <  Cklt1  -  tQ|  1/,J  for  each  k  and  tQ, tt  0  (p.T^. 

Lemma  3.1:  There  exist  bounded,  continuous  functions  (vn(x, t) ,wn(x,t) )  which  satisfy  the 
Equations  (3.1)  in  ST^. 
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Proof:  The  proof  follows  Evans  and  Sheni:  [1].  Let  v_  (x,t)  -</>(x)  and  w_  (x,t) 
-  *0  n0 

w(x,0)  =  0  in  Sm  .  Assuming  that  v„  (x,t)  and  w„  (x,t)  have  been  defined  for 
M  nj  nj 

j  >  0  ,  we  let  vn  <x,0)  -  ^(x)  and  w_  <x,0)  •  w(x,0),  and,  for  (x,t)  6  Sm  , 
nj+1  nj+1  A1 


(3,2a) 


os  t  n  ' 

v  (x,t)  -  /  K(jc-5,t)v><e)dc  +  r  at  /  K(x-s,t-T)dc 


0  -s  (t) 
n 


-  Jfc  dT  f  K(3t-e,t-r)w  (c,T)de 


0  -» 


w  (x , t )  -  ee  eYt  J*  eETTV  (x,nMn 
Vi  0  nj 


The  resulting  sequences  of  functions,  {v  <x,t) }  and  {w  (x,t)},  are  defined  and 

"j  nj 

continuous  in  ST^ •  We  show  that  these  sequences  converge  uniformly  to  a  solution  of  the 
Equations  (3.1).  Note  that  since  1(x)  is  bounded,  it  follows  from  induction  that  each  of 
the  functions  v„  and  w„  are  bounded « 


p,  (t)  -  sup  { I v  ( x, t )  -  v  (x,t)|  4*  | w  ( x, t )  -  w  (x,t)|). 

3  (x,t>es  nj  nj-i  "j  nj-i 

1 

Fro*  Equations  (3*2)  it  follows  that,  for  (x,t)  e  , 

( 3*  3a)  I  v  ( x,  t )  -  v  ( x ,  t )  |  <  /*  dx  f°  K(x-£,t-i)|w  (£,t)  -w  (£,T)|d£ 


(3.3b)  I w  (x,t)  -  w  ( x , t )  I  <  e  J*  |v  (x,n)  -  v  (x,n)|dn  • 

Vi  nj  o  nj  Vi 

Adding  (3.3a)  and  (3. 3b)  we  obtain 

p.  +  1(t)  <  (1  +  e)  J*  p ,  ( T ) d t  for  t  e  <0fTt),  j  -  1,2, ••• 
J  0  ^ 

If  K_  is  a  bound  on  p^(t)  for  0  <  t  <  T^  we  have 


Pl(t)  <K  ,  p2(t)  <  (1  +  c)Kt,  < 


K[(1  +  e)t] 


jii  V  -ji,  <  j+i  >  i 


...  .J+l  ( 1+e)T 

1 : ,  -  <  Ks  and  the  sequences  {v  (x,t)}  and 

( j+1 ) I 


{w  (x,t)}  converge  uniformly  in  S-  to  limit  functions  vn(x,t)  and  w  (x,t). 

1 y  n  n 


Moreover,  passing  to  the  limit  in  equations  (3.2)  we  find  that  vn(x,t)  and  wn(x,t) 


satisfy  the  Equations  (3.1)  in  ST^.  // 

Prom  the  proof  of  the  preceding  Lemma  it  follows  that  there  exist  constants  V  and 

W  such  that  |v  (x,t) |  <  V  and  |w„(x,t) |  <  W  in  S«  .  Note  that  V  and  w  can  be 
n  n  i-  j 

chosen  Independent  of  the  curve  sn(t).  From  (3.1)  it  follows  that  ^vn^x^x#t^  and 
(wn)x(x,t)  both  exist  in  ,  except  possibly  at  |x|  -  sn(t),  0  <  t  <  . 


Lemma  3.2.  There  exist  constants  and  ,  independent  of  the  curve  Bn(t)  ,  such 

that  |(vn)x<x,t)|  <  V1  and  |(wn)x(x,t)|  <  W1  in  ST^,  except  possibly  at  |x|  -  sn(t), 

Proof:  Suppose  that  x  *  s_(t).  We  differentiate  both  sides  of  (3.1a)  to  obtain 
-  n 

_  *  »_(T) 

(vn)x(x,t)  -  /  Kx<x-5,tK>U)dS  +  r  dT  /  Kx(x-e.t-T)dC 

-«  0  -s  ( T) 

n 


-  /fcd T  r  Kx(X-C,t-T)Wn(C,T)de  . 

Integrating  by  parts  in  the  first  integral  yields 

Kvn)x(x,t)i  <  iv>,(*)ia>+  (i  +  wj  /fcdT  /”  iKx(x-e,t-T)ide  . 

*  0  —« 

Note  that 


^dx  /*  |K  (x-C#t-T)  |d£  -  2  /*  dt  ( Tift-x)  |dn 

0  -oo  o  0  * 


-(t-T) 


7% 


t-T) 


dT  f 
0 


2(  t-x) 


e’n  /4(t"T)  dn 


J*  TX - Ti  dT 

/2(t-x)  /2 


0  w 


Therefore, 


<  T 


1 


v2. 


|<v  )  (x,t)|  <  l/(x)  I 
n  x  * 


+  (1  +  W)T4 
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.  -  * 


Differentiating  both  aide*  of  (3.1b)  yield* 

l<»n)je(**t)  I  <  c  J*  l(vn)x(x,n)  Idn 


<  c  V1  =  W1  .  // 


ia  *  bounded  continuous  function  in  8^ 


Proof :  follow*  because  (w)  -  c(v  -  yw  )  in  8*  .  We  choose  W„ 

nr  n  n  2 

l(Vt(*,t)|  <  W2  in  8,  .  // 


Let  W  -  *»,  ♦  \  *  »3  • 


(3.4) 


W.  wish  to  dafin.  •n*.j(t)  Implicitly  as  th.  solution  of  th.  aquation t 
*(,n+i(t)'t)  "  •  "  *<•  Ht)  +  r<w  )<s  (t),t)  , 


W0)  "  x0  • 


Rscall  that  wo  aro  assuming  that  sn(t)  is  a  smooth  function  in  (O.T,), 

M  -  sup  s‘(t)  <  •  ,  and  given  p  e  (0,T,),  th.ra  exists  a  constant  C 
t«(0,T.)  n  1 


i*;<v  ■  4  cn,trtol  ^  for  ,aoh  *0**1  ®  <p>v  prom 

we  conclude  the  following. 

bat  8(t)  equal  the  right  hand  side  of  (3.4).  Then, 
a)  8(t)  C1(0,T1)  , 


2.3, 


b)  there  exists  constants  Ky  and  K0  such  that 


(3.5) 


I  B1  ( t)  I  <  k7  ♦  in  (OpT^  , 


c)  there  exist  constants  Kg  and  K10  such  that 

le’(V  -  B*<V'  <  (K9  +  K10  Cn  T11/2)  1/2 

for  each  ^  ,tQ  0  (p, T^. 

Note  that  the  constants  Ky  and  Kg  depend  only  on  W  ,  and  are,  therefore t 
of  n  •  Furthermore,  Kg  and  K^q  depend  on  p  and  the  bound  on  |0'(t)| 


so  that 


such  that 

2.4  and  2.5 


independent 
given  in 


(3.5b).  Hence,  Kg  and  K1(J 


can  be  chosen  independently  of  n  • 


We  conclude  from  Lemma  2, 1  and  the  Implicit  function  theorem  that  there  exists  a 
smooth  function  sn+1(t),  defined  for  some  time,  (say  t  e  [0,T21),  which  is  a  solution  of 
(3.4).  We  show  that  as  long  as  (sn+1(t)#t)  stays  in  the  rectangle  P,  defined  in  Lemma 


2.1,  then  s^+^(t)  is  bounded,  independently  of  n 
We  differentiate  Equation  3.4  to  obtain 


Vn+i<t>'t),;*i(t)  +  -  e’(t)  - 


(3.6) 


From  Lemma  2.1  and  (3.5b)  it  follows  that  if  (sn+1(t),t)  6  P,  then 

<T  [K7  +  K8MnT1/2+  S3] 


K11  +K12  MnT1 


v2 


where  ^  -  —  (K^  +  5^)  and  -  1^/6^  do  not  depend  on  n  . 


12  8  2 


Suppose  that  <  (  — 1  )2  .  Then,  as  long  as  C ^  ( t: >  # t )  e  P, 
**12 


Hence , 


|s;>i(t)|  <Kn  +1/2m„  • 

Kn+1  «  *11  +  V2Mn  «  *11  +1/2‘K11  +  1/2M„-1>  «  ***  4 


<  *,,<»  +  1£  +  •**  ♦  >  ♦  2""  M0  <  2*,,  +  M0  =  H  . 

Therefore,  the  sequence  {s^(t)}  is  uniformly  bounded  by  the  constant  M  .  It  follows 

that  there  exists  a  constant  T  such  that  T  <  T ^ ,  and  (sn(t),t)  e  P  for  each  t  e 

(0,T)  and  each  n  .  Furthermore,  there  exists  a  subsequence  {s  (t)}  which  converges 

nj 

uniformly  on  (0,T]  to  a  continuous  function  s(t).  We  assume,  without  lose  of 


generality,  that  (s  (t)}  -  {s  (t)}. 

n  j  n 


such  that 


Lemma  3*4:  Pix  p  6  (0,T).  Thera  exist  positive  constants  and 


n+ 


1**1*  "  *  *  (K13  +  *14  Cn  T1^  *  *fc1  ”  2 


for  each  n  and  t^t^  9  <p,T)  .  The  constants  and  can  be  chosen 

independently  of  n • 

Proof :  This  follows  from  (3.5c) ,  (3.6),  and  Lemma  2.1.  // 


We  now  assume  that  T  <  (  ^ •  Then  the  previous  lemma  implies  that 


14 


<  +  V2C  <  •••  <  K  ,<1  +  %*  •••  +  — TT)  + 


2  H. 


That  is,  given  p  9  (0,T),  there  exists  a  constant  C  such  that 

-  »i<vi  <  c  it1  -  tQi 1/2 

for  each  n  and  tg,tj  6  (p,T).  It  follows  that  s*(t)  is  continuously  differentiable 

in  (Q,T)  and  a  subsequence  of  (sMt)}  converges  uniformly  on  compact  subsets  of 

( 0 ,T)  to  s'(t).  With  loss  of  generality  we  assume  that  {s'(t)}  converges  uniformly  on 

n 

compact  subsets  of  (0,T)  to  s'(t).  // 


Lemma  3.5*  The  sequences  {v  }  and  {w  }  converge  uniformly  in  to  continuous 

-  n  n  * 

functions  v  and  w  which  satisfy  the  Equations  (1.4). 

Proof:  bet  p  (t)  -  sup  {|v  (x,t)  -  v  (x,t)|  +  |w  (x,t)  -  w  (x,t)|)  .  Prom  (3.1a) 

-  n  n+ 1  n  n+ 1  n 

705WL 

it  follows  that  for  (x,t)  0  ST  , 
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I 


*  8  i  i  (  t) 

|vn+1(x,t)  -  vn(x,t)|  <  I  /  dr  /  K(x-£»t-T)d£  - 


/Sx  /  "  K(x-e,t-T)de|  +  /Sx  /“  K(x-C,t-X)  |w  ( C,  T)  -wn(e.x)|de 

0  -s  (x)  0  — 

n 


*  "8  (t)  Jr  ,8n+1<T) 

|  j  dx  /  K(x-e,t-x)de  *  r  <i x  /  K(x-e,t-x)dd 


0  -Vi(T) 


0  *  <T) 

n 


+  j^dx  f°  K(x-5,t-x)  |wr+1  (  C#  t)  -  wn<C#T)|dC 

0  —  OD 


<2  sup  Is  (x)  -  s  <x)  I  J* - 17~ - l/,dT  + 

0<x<t  *H  "  0  2x/2<t-x)/2 


+  J*  sup  |w  (x,  x)  -  w  (x,x)|dx 


(3.7)  <  2T  ^2  SUp  |s  (x)-s  (x)|  +  /Sup  |w  ,(x,x)  -  w  (x,x)|dx  . 

_  n+  in  i  __  nt  i  n 


Prom  (3* 1b)  it  follows  that  for  (x,t)  e  0T  , 

(3.8)  |w^+1(x,t)  -  wfl(x,t)|  <  c  /t|vnf1(x,T)  -  vn(xfx)|dx 


Let  6  -  2T  '2  sup  |s  .  (x)  -  *  (t)|.  Note  that  6  +0  as  n  ♦  «*.  Adding  (3.7)  and 

n  n+1  n  n 

0<t<T 


(3.8)  we  find  that  for  t  6  (0,T), 


p  (t)  <  6  +  (1  +  e)  /t  p  (x)dx 


,r^ ^ m  ’ 


I 

I 


Proa  Gronw.ll*.  Inequality  It  follows  that 

p  (t)  <  C  « 
n  n 

for  iom  constant  C  independent  of  n  .  Therefore  p  (t)  ♦  0  uniformly  as  n  ♦  •  ,  and 

n 

the  sequences  {v  }  and  {w  }  converge  uniformly  to  limit  functions  v(x,t)  and 
n  n 

w(x,t).  Passing  to  the  limits  in  <3.1)  we  find  that  ( v(x,t) ,w(x,t) )  satisfies  the 
Equations  (1.4).  This  implies  that  v  and  w  are  continuous  functions  in  ST  •  // 

Theorem  3.6:  Let  K  -V4«in<  V2  -  a, a)  and  suppose  that  T  <  — ^  .  Then  (v,w)  is  a 
classical  solution  of  the  Cauchy  problem  (1.1)  in  ST  • 

Proof ;  Throughout  this  proof  we  assume  that  t  0  (0,T).  Recall  that  v( x,t) ,w( x,t ) ,  and 
s(t)  satisfy  the  Equations  (1.4)  and  (1.6).  Setting  x-e(t)  in  (1.4)  and  subtracting 
the  resulting  equation  frosi  (1.6)  we  find  that  v(s(t),t)  -  a  • 

Equation  (1.4)  implies  that  for  |x|  <  s(t),  (v,w)  satisfies  the  differential 
equations 

v  -  v  -  v  +  1  -  w 
t  xx 

(3.9) 

wt  -  e(v  -  yw) 

and,  for  |x|  >  s(t),  (v,w)  satisfies  the  differential  equations 

\  *  Vxx  "  V  "  * 

(3.10) 

w  -  e(v  -  yw)  . 

We  show  that  v  >  a  for  |x|  <  s(t),  and  v  <  a  for  |x|  >  s(t).  This  implies  that 
for  x  *  s(t),  (v,w)  satisfies  the  system  of  Equations  (1.1). 

Suppose  it  were  not  true  that  v  >  a  for  |x|  <  s(t),  and  v  <  a  for  |x|  >  s(t). 
For  example,  suppose  that  v(Xj,tj)  <  a  where  |x^|  <  s(t^).  Since  v(x,0)  >  a  for 
| x |  <  Xq  -  s(0),  we  may  assume  that  v(x1,t1)  -  a  and  v(xft)  >  a  in  the  region 
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\ 


G  -  {{x, t )  1 1 x |  <  a(t)#  t  e  (0,^)}.  We  use  the  maximum  principle  (see  [6),  page  159)  to 
show  that  this  is  impossible.  Note  that  v  -  a  for  |x|  -  s(t)  and  v(x,0)  >  a  for 
|x|  <  Xq  •  Let  L  be  the  operator  defined  by  I 2  v  -  +  v  •  Then,  in  G  , 

Lv  -  1  -  w.  From  (1.4b)  it  follows  that  in  R  x  (0,^)  , 

|w(x,t)  I  <  e  |v(x,n)  Idn 


Therefore,  in  G,  Ur  >  1  -  K  >  a  -  L(a).  It  now  follows  from  the  maximum  principle 
that  v(x1,t1)  >  a  .  This  is  a  contradiction.  A  similar  argument  shows  that  it  is 
impossible  for  v  >  a  for  | x |  >  s(t). 

We  have  shown  that  except  for  x  *  s(t)  ,  (v,w)  satisfies  the  system  of  Equations 
(1.1)  in  ST.  It  remains  to  show  that  vx(x,T)  exists  for  |x|  -  s( t) • 

Assume  that  |x|  <  s(t)  and  |£|  <  s(t).  Then  ( v( £, x) ,w( x) )  satisfies  the  system 
of  equations 

v  -  v  +  v  -  1  -  w 
T  ££ 

w^  *  e( v  -  yw)  . 

Multiply  both  sides  of  the  first  equation  by  K(x-£,t-x),  integrate  by  parts,  and  use  the 

fact  that  K  +  K  _  -  K  -  0  to  obtain i 

t  CC 

(Kv)t  -  <Kv^  +  (K^V^)  “  (1  "  W)K  • 

We  integrate  this  last  equation  for  -s(x)  <  £  <  s(x),  6  <  x  <  t  -  6,  and  let  6  +  0  to 


obtain: 


TT3 


X  - 

v(x,t>  -  /  X(x-$,t)*U)d£  -  /  K(x-»(T),t-T)a*'(T)dT 


-  J*  K(*f»(T)  (T)dT  -  /*  K(X-*(t)#t-T)V  («(T)’rT)dT 


(3.12a) 


+  K(x+s(t)  ,t-t)v  (-s(t)  +  ,  T)dT  ♦  a  K  (x-i(t)  ,t-T)dT 
0  5  0  * 


-  /*  a  K  (x+a(  t)  ,t- x)dx 


h* 

% 


i 


/*  «jt  /s(t)  (i  -  %*)K<*-e,t-T)d5 
0  -s(  r) 


Next  assume  that  £  >  s(t).  Then  v(£,r)  satisfies  the  differential  equation 

Vt  “  VCC  *  V  “  # 

Multiply  both  sides  of  this  equation  by  K(x-£,t-T)  and  integrate  by  parts  to  obtain: 

(Kv)  ^  -  (Kv^^  +  (Kev)c  *  “Kvr  " 

Integrate  this  last  equation  for  s(t)  <  £  <  »  ,  6  <  t  <  t  *-  6  and  let  6+0  to 

obtain: 

K(x-£,t)*(£)d£  ♦  /fc  *(*-•(  T),t-T>*»*  (T)dt 
(3.12b)  +  K(x-.(t)  ,t-T)v  (s(t)  +  ,t)<1t  -  /C  a  K.(*-s(T)  ,t-T)di 

ft  *  ft  * 


Jfc  dT  J*  K(x-C.t-T)w(  5,T)d5 
0  s(t) 


Similarly,  for  5  <  s( x)  we  obtain: 


/  °  K(x-C#t-T)dC  +  /*  K(xfS(  T),t-T)a  S*  (T)dT 

—os  0 


(3.12c)  -  J1  K(X+s(t)  ,t-T)v  (-s(  T)  ,T)dT 

ft  ** 


+  J*  a  K  (Xfs(T)  ,t-T)dT  -  l^dT  f“fl  T)K(x-C,t-T)w(CfT)dC 


0 
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Adding  {3.12a),  (3.12b),  and  (3.12c),  and  using  (1.6)  we  find  that 


(3.13)  T)  ,t-x)  (Vr(s(T)+,  T)  -  Vr(*<T)",T)J 

0  *  * 

+  K(XfS(T),t-T)(v^(-S(T)+,T)  -  V^(-s(x)  ,T)])dT  •  0  . 

Using  the  assumption  that  y?(x)  -  <^(-x)  it  follows  from  (1.4)  that  v(x,t)  -  v(-x,t)  in 
ST  .  Therefore,  (3.13)  can  be  rewritten  as 

/t(K(X-S(  T),t-T)  -  K(x+S(  T),t-x)]  (vf(s(  x)*,x)  -  V  (s(x)  ,T)]dt  -  0  . 

0  *  4 

Since  [K(x-s( t)  ,t-x)  -  K(x*s( x)  ,t-x)3  >0  in  (0,T)  we  conclude  that  vx(s(t)“,t)  - 
vx(s(t)+,t)  in  ( 0, T) . 

me  have  shown  that  vx(x,t)  is  a  bounded  continuous  function  in  ST*  Prom  (1.4b)  it 
follows  that  wx(x,t)  is  also  a  bounded  continuous  function  in  ST  •  // 


-2  1- 


REFERENCES 


Evans,  J.  and  N.  Shank,  Solutions  to  axon  equations,  Blophys.  J.  10,  (1970), 
1090-1101. 

FitzHugh,  R. ,  Impulses  and  physiological  states  in  models  of  nerve  membrane, 
Blophys.  J.  1,  (1961),  445-466. 

Friedman,  A.,  Partial  Differential  Equations  of  Parabolic  Type,  Prentice  Hall, 
Englewood  Cliffs,  N.  J.,  1964. 

McKean,  H.  P. ,  Nagumo's  equation.  Advances  in  Math.  4,  (1970),  209-223. 
Nagutno,  J.,  S.  Arraoto  and  S.  Yoshizawa,  An  active  pulse  transmission  line 
simulating  nerve  axons.  Proc.  Inst.  Radio.  Eng.  50,  (1962),  2061-2070. 

Protter,  M.  H.  and  H.  F.  Weinberger,  Maximum  Principles  In  Differential 
Equations,  Prentice-Hall,  Englewood  Cliffs,  N.J.  (1967). 

Rlncel,  J.,  Impulse  propagation  in  excitable  systems,  in  Dynamics  and  Modelling 
of  Reactive  Systems  W.  K.  Stewart,  W.  H.  Ray  and  C.  C.  Conley,  eds.) ,  Academic 
Press,  New  York,  (1980),  259-289. 


/db 


iCCjRlTY  CLASSIFICATION  OF  THIS  PACE  fl»7u*n  fnteted} 


j  REPORT  DOCUMENTATION  PAGE 

K K A l>  INSTRUC  TIONS 

MF.KORK  COMPl.f  i  INC,  FORM 

1.  REPORT  NUMBER 

2221 

4.  TITLE  (end  Subtitle) 

LOCAL  EXISTENCE  FOR  THE  CAUCHY  PROBLEM  OF  A 
REACTION-DIFFUSION  SYSTEM  WITH  DISCONTINUOUS 
NONLINEARITY 

5.  TYPE  OFJtEPO^T-ft^CRIOO  COVERED 

Summary  Report,  -  no  specific 
reporting  period 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTNORp; 

David  Terman 

•  // 

8.  CONTRACT  OR  GRANT  NUMBERf#; 

DAAG29-8O-C-0041  ’ 

MCS  8  0-17158 

9.  performing  ORGANIZATION  name  and  address 

Mathematics  Research  Center,  University  of 

610  Walnut  Street  Wisconsin 

Madison,  Wisconsin  53706 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

1  -  Applied  Analysis 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

See  Item  18 

*2.  REPORT  DATE 

May  1981 

13.  NUMBER  OF  PAGES 

22 

1 4,  MONITORING  AGENCY  NAME  4  ADDRESSfi/  different  from  Controlling  Office) 

1 

IS.  SECURITY  CLASS,  (ol  (hi,  report; 

UNCLASSIFIED 

1S«.  DECL  assi  fi  cation/  downgrading 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  (hi,  Report) 

Approved  for  public  release;  distribution  unlimited. 

t7.  DISTRIBUTION  ST  ATEMENT  (o t  the  ebetract  entered  In  Block  20,  it  different  from  Report) 

1®.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office  National  Science  Foundation 

p.  0.  Box  12211  Washington,  D.  C„  20550 

Research  Triangle  Park 

North  Carolina  27709 

19.  KEY  WORDS  (Confinue  on  revecae  side  if  neceaaary  and  Identify  by  btock  number)  i 

Reaction  Diffusion  Equations,  Classical  Solution 

0 

?0.  A9STRACT  (Continue  on  nreree  aide  it  neceeeary  end  identity  by  block  ntiiiber) 


The  pure  initial  value  problem  for  the  system  of  equations 

vt  =  vxx  +  f(v>  “  w 
w  »  e (v  -  yw) 

is  considered.  Here  e  and  y  are  positive  constants,  and 


DD  ,'anM7J  1473  COITION  or  INDV«  IS  OBSOLETE  UNCLASSIFIED 

SCCiTriTY  CL  ASs7fiCAT|ON  OF  THIS  PAGE  (WTten  ~l~nt*r*d) 


20.  ABSTRACT  -  Cont'd. 

f (v)  -  v  -  H(v  -  a)  where  H  is  the  Heaviside  step  function  and 
a  e  (0 )  •  Because  of  the  discontinuity  in  f  one  cannot  expect  the 
solution  of  this  system  to  be  very  smooth.  Sufficient  conditions  on  the 
initial  data  are  given  which  guarantee  the  existence  of  a  classical  solution 
in  1R  *  (0,T)  for  some  positive  time  T. 


